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PREFACE 


The purpose of this contract was to provide data processing 
and analysis in support of two 6SFC programs: tornado detection by 
analysis of radio frequency interference in various frequency bands 
and sea state determination from short pulse radar measurements. 
Extensive data reduction efforts were performed in support of these 
programs. A new and unique backscatter simulation was implemented 
to predict radar performance as a function of the wind velocity. 
Computer programs were developed for the various data processing 
and analysis goals of the contract effort and delivered to 6SFC. 



I. Introduction 


The data supplied to Techno«Sciences for reduction and analysis 
consist of 9 track, 800 bpi digital tapes in various record formats 
and include, in addition to the data records, records indicating the 
wind state at the measuring platform, aircraft orientation, altitude, 
and time of recording. These data were collected over several 
experiments and have several corresponding formats. Therefore it was 
necessary for Techno-Sciences to develop highly flexible programs to 
be able to handle the various formats. Over 60 tapes have been involved 
in the analysis and data reduction job. These tapes were generated 
under a variety of sea conditions providing an effective data base 
for investigation. 

An important tool in the difficult task of determining sea state 
conditions from short pulse radar data is the use of simulated data. 

Data generated by simulations is highly controlled so that it is 
possible to find and evaluate estimation methods without uncertainty 
as to the true value of the quantity to be estimated. Techno-Sciences 
devloped a 2 dimensional computer simulation during the contract period. 
The modelled data has a high degree of correlation with the measured 
data in appearance. 



II. Short Pulse Radar Simulation 


The problem of sea state estimation from measured short pulse 
radar information is very complex. Particularly complicating the 
problem is the lack of accurate "ground truth" data, i.e. the lack of 
exact knowledge of what an accurate estimate on any particular segment 
of data is. A good way of removing this problem is to develop computer 
simulations where the estimated parameters can be accurately known and 
controlled. 

As a part of the contract effort, Techno-Sciences has developed 
a two-dimensional simulation of the received pulse under a Gaussian 
surface model. A listing of the program appears in Appendix A along 
with a sample printout of several pulses. The surface is allowed to 
consist of any combination of a deterministic sinusoidal component 
and a Gaussian component with a Pierce-Moskowitz spectrum under operator 
control. In each case, the period and amplitude is variable. The random 
portion is generated using an FFT on randomly generated frequency com- 
ponents whose variancsare porportional to the Pierce-Moskowitz spectrum 
at each frequency. The simulation of the pulse return is based on equa- 
tion (4) of Spectrum of Power Scattered by a Short Pulse From a Stoch- 
astic Surface , by David Levine, NASA report X-952-74-299, August 1974. 
The specular reflection points are found by interpolation between sam- 
ples of the derivative of the surface where it goes through zero. 

It has been found that the simulated results agree well with the 


measured data. 



Appendix A 

SHORT PULSE RADAR SIMULATIOH 
Program listing and Sample Printouts 



Program SIMUL - Program Listing 

fASSM 

SIMUL FRSd 2 D SH'^FT PULSE SLhTTER SIMULA! I ( 5 N LDD 
EFL5T 
fPORT 

C FR0GFAM T0 SIMULATE A liOMrlCATItjN !2 jF EON (4) IN THE TIME 

c dumhin rr •'Spectfum &f power — ■■ &r levine 0F aug 1974 

C THE sThTIONHFY POINTS ArE F0UNP Fk)R A GIVEN SURFACE 
C USING FFT TECHNICUES 

C EOTH A SINUSdlDAL hND A GhUSSIAN SURFACE 
C COMrDNENT ARE All^V ED T nE SukFhCE 
C IS GENEPHTEr IN THE FREOUENCt D0MH1N 
C TO r^ET THE DESIRED SPECTRUM THE 
C RhNDuM SURFhCE IS FIEF' E-h'jSK0WITS 
C WITH GIVEN PEAK FPEQUEN.Y AfJD 
C RMS ' KLUE as INPUT THE SINUS&IDkL 
C FhkT hhS THE ^.hriE PERIOD buT ANY 
C 0THEF CHOSEN RMMS VALUE 
C FR'jial hM by TECHNai-:CIENCE5.INC 
r 

C 

IMPLICIT INTEGER+2 rl-N^ 

INTEGc.R4>^ iRkN 

COMFlE> 1DF4 ',5m 1024 g SI'DC 102-^ « CMPLX, C0NJG 
COMFLEa W, 5S 

DIMENSION TEU. S 13 IgFD' El ' lFLnT(2n LASTtSlJ 
DIMENSION FhrGcf 1024 ^-ODT^ .U24 •) 
t-CUi '/ ALE^^CE ' FhN'IE' 1 1 ' ' ' CR''/^'C'T^. I t 2DIh I 

T hrK < >=SIN* A '•r rSf ^ ) 

C SFE:( ) = FIERSON-rM^ h ij .its SFECTPLm 

EFE*>a)=H 24E =(-^MlN.riro . E G2U4/- C X* X ) '^< + ♦3'' 

C ^hNGliLHF EJFFhCE FFE'OUENrY 
C IFnN=FAfiD0M NlfiEEF 'IENE-hT'jF 
C GhUSdIhN PULSE zFhFE 

FULSE( /^=E.'Ff“W '^.'S.IGFLS 3 

C ' SItk 3 ^ ’♦♦2 =hNTE'(Nh FUV*E=ChIN FATTEFN 

GrIN' V j=. SINi ^^H.jG‘/rtL j.hMh/K A3 Sk \,"ANG'^hL ), 1 E-S) ^*»K2 
ThTh SIGTIM Quo/' 

ZATh IRAr</'lc7 1,2l,'' 

DhT a NF”T /• i C 2 4 
L'ATH LU^S'' 

DhTA P2 =' 2:r .53 ' 

data C,DT FULSnI. FO/300 00-25 ,01, 13500 z' 

: C^EFEED OF LIGHT DT = TIME I i' F , FIJLE TM= F UL5E DURATION 
IATA I3/3> 

Data 14/4 

C DTuR =*'EGRElS TCi FhDIhNS C '^uD' L^ 51 0 J 
rPTH Dior/ Ol^‘i-323'' 

DhTh I0NE,IEIg, 1, 30000/ 

I-hTh IGDriS^'x 
LC'=F0*F2 
r.FFT2=NFFT 2 
f^FFT 2 1 = NFF7 2 1 
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WRITEiO 25) 

25 Fi£iPriAT( 53HENTEF RM5 SINUSOIDAL'' M ). 

I RMS RANDidMCF-n*ai PERI0DC M )fF8) } | 

FEAD(0, 2&; PMSSIN, PMSRAN, FEPI0D , 

2^ F0RMATC SF8 > 

WRIT EC 0. 27 ) 

27 FORMATC 33HENTEF RC 3DB CUTOFF^PULSE DUR( FS ) ) I 

REHDc 0- Z6 ) FCTS 
UFITEC0,31 ) 

31 FORMAU 30HWAVE REPEAT S. PULSE REPEATSCIS)) 

PEAlKO. 32; IRPTW, IRPTP 

32 FORMhTC 815 ) 

WRITE' 0. 41 ) 

41 FORMAK llHBLlZiCI<S'?CI5;) 

FEADcO-32; IcLkS 

FH0=EXP( -P2*F0TS*DT/FULSTM; 

PhOUSORTt 1 -RHB*RHO ; 

BG2LI4-F2^P2* 1 5/f PERI PERI DD ) 

A2 = P2»RMSRArURMSRAN4iEG2U4 
ALTUD=2000. 

C CON-'ERT FROM PULSE TIME SIGMA TO RANGE SIGMA 
SIGPLS*SIGTIM*C2 

C RNGINC^ASSUMED PULSE STAFT.END POINT 
FNGINl*^ *SIGPLS 
ANGLIN* 18 

ANGCTR*ANGLir.iDT0R 

BMWITH*2^ 

HNGIDC=5r';WITh4DT0R/2 
HNGMh \= AN‘jCT F f f-.r<G3Db 
hf.'Gni N= hNGr IF- HNG3D5 

(-rJGVAL*HNGIDE 1 33 I 

YD=HLTUD4TANC HNGMIN ) I 

C SET L'F GRID 

GDINC=13SS .FLi^iAU IGDFT5- 1 J ! 

DU 5 I = 1 * IGDF T S I 

E IGFDc I ;*FLOhTi I- 1 ;*GDINC +.5 | 

Y3DB=ALTUD4( TANC Ar^iMHX ^-TANC AN':,MIN ^ ) 
pr=irrDB rL-jnT' nfft ; 

DY2 = Dt^-2 
Y02=tQ-DY2 
0 DF+Dt= l^NFFT 

DF= 1 /Y3DB 
FTDF = SuRT( DF . 

DW=P24DF 

C D'j=DELTh uMEGh, DF = DELTH FREQ 
C SET UP FFT TAELE 

call FOUR IT' 5 NFFT, TEL, O; 

C SET BLOLV C0UIJTEF T0 EERO BEFORE L'-^OPING 
IC0UNT=0 

ET UF RhNDOM fart OF WAVE TPANSFUFM WITH SFECTPUM SFECcX) 

C FhRT5=0 
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40 2t 1 )^C\1FLXk 0 , Cl ) 

2D( 1 )«CMPLA( 0 > 0. ; 

2DP< i )=CriFL^< 0.0.) 

D0 1 NFFT2I 

X = FL0AT( I- 1 )*D\^ 

H( I )=CMPLX( BORIC 5FEC( X ) IRAN ), 0. )*RTDF 

1 CONTINUE 
IFREQ«T3DE/P^EKI(3D + . 5 

C SET UP N 0 N-RANDCri WAVE C0MP0NENT 

2t IFREO )=2C IFREQ )-*>CMPLX^ RMSBIN/^QRTC 2 ). 0. ) 

C SET UP F0R PL0T SCALE FACTORS 
SCHLE»RnS5IN**2 
X = FL0ATC IFREO- 1 ;*DW 
SCALE 1=SCALE*X»X 
SCALE2=5CHLEUX*X 
1)0 2 1-2. NFFT21 
X = FL0ATc I- 1 ;»DVf 
XX = SFECC X ♦I'F 
5CALE«SCALE+>X 
XX=X<#X*X 
SCALEl=5CALEl4>X 
^X«XX*X*X 
SCALE2=3CHLE24>X 
W=ChPLXC 0 . X .> 

HDt I )«W*EC I ) 

EI^Dk I )=W>HD( I ) 

2 C0UTINUE 

5crle=^bqrt( scale > 

BCHLE1 = SGFT( SCALEl ) 

SCHLE2*5QFT>' SCHLE2 ; 
h\S 3 I = 2.NFFT21 
lL=NFFTf2-I 
Ef IC )«CONJG( Ei' I ) ) 

= D( IC )=CONJG< EIU I ; ? 

EDDc IC '>=C10NJGC EDDC I ) ) 

3 C'JNTINUE 

CALL FiOURlTv E, NFFT, TEL I0NE ; 

CALL FOURiU ED. NFFT TEL. lONE) 

CALL F0UP m EDI'. NFFT. TEL. lONE ) 

C FIND STATIi^tNHFT F0INT5 ASSUMING LINEAR INTERPOLATION 
C hfiD PLOT Hr r ' FIF5Tc_.2D DE^ I ATI E 5 

Call SETGrD). lU. iGr D I'jDFTS ItrIG ) 

TMDIF=C Y02 + Dt ALTUD - fehlcedci;; 

I” sO 

DO 4 1=1 NFFT 
t«FL'>iATt I >*Dt + T02 
THNTHT= Y/ALTLD 
TPNhLF=REAU EDc I ) ) 

TANDIF=TANTmT-T analf 

IF C THNDIF*TriDIF. GT. 0 ) G0 TO ^ 

DLT = TMDIF^^ TMDIF-THfJDIF ) 

DLT1=1 -DLT 
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DLT*DY 

HEIGHT*FEhL( r( 1-.1 ) )*DLT UDLTtFsEHU 2( I ; ) 

DERIVoREriLC 2Dc I- 1 ) )*DLT l4 DLT ♦ T qNr^LF 
DERIV2=RErtL< 2DI'( I - I ) ;*DLT l^IiLT^REPlLC 2DDf I ) J 
K«K+ 1 

RfiNGEt. K )=S0RT( ALTUD-HEIGHT )**2 ) 

HNG«ATflN( PERIV • 

FADCRV«U 1 4lEFIV*DERIV 1** 1 5 )/PERIV2 

CR00T^ K )eSQRT( hBSC RhDCR^A ^ RHDCRV-RANGEi K ) ;*RANGE^ K ) ) J ) 
CR00K K )=CR00TC K )»GAINC ANG-ANGCTR > 

6 C0NTINUE 

IF( IRPTU. EQ 0 )Ci0 T0 4 

IFL0TC 1 )«20. *RERL( 2( I ) ;/5CALE + 233. 5 

IFL0TC 2 )»20. IRERLC 2DC I ^ )^5CRLE 1 + 700. 

IPL0T(4) = 2O ♦TRNTHT.'5CALE1 + 7D0. 

IFL0T( 3 ;«=20. *RERU HDf( I ) J/SCHLE2+ 1 U7. 5 
D0 30 J=1,IRPTV/ 

C REPEAT T0 SPREAD SCALE 

30 CALL PL0TITCIPL0T, LAST. 14, I0NE ) 

4 TliDIF-TANDIF 
NUM»K 

C FE0RDER IN INC«^EASING RANGE 

10 DO 7 K»2. NUM 

IF' RANGED K- 1 ) OT.RANGECK)) GO TO 6 

7 CONTINUE 
G0 T0 9 

5 X=RANGE' y ^ 

FhiNGEt K )=Fhr GEC k- 1 > 

FhrJGEC K- 1 

LPO0T < K '' 

CFOuTr k i^CFuC'^f K- 1 ;• 
r'FUOTC k- 1 '=A 
GO T0 10 
,0NTINUE 

GEt.E-nTE 0UTFUTE NUi/ 

TItlE=( RAN>:,Ei 1 j-RNGIN: >/'C2 
Y*0 

I hST« 1 ^=0 
IFLO'f 1 ;=o 
u'' I T E*" 7 L 1 ) 

2 1 FCFflHT ( ImO > 

11 TinE=TinE^DT 
F G = C2* T IMl 

C RG=RANGE F0R FULSE CENTER 
c rgmh ,rgmin=sthft end 
PGMAX=RG+PNGINC 
FGhlN>=Rb-RNGIN( 

C FIND FEFLECTOR Liril TS. . MAX'' 

C I E Thki5E IN THE F hNGE BIN 
DO 12 1=1. NUN 

IFcRhNGE(I) GE RGniN; GD T0 i3 

12 CONTINUE 
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ICOUNT=IC0UNT^ 1 

IF(IC0UNT LT lELKS) GO TO 40 

STOP 


8 

13 

f1INV*I 

DO 14 lorilNV-NUM 

IF( RANGE( I GT. RGliRO G0 TO .5 



14 

CONTINUE 
Ii^NUIV 1 



15 

MAXV = I- 1 
2H=CMPLX( 0. . 1 



C GENERATE 0UTPUT3 
C iL PLOT THEM 




IFCMAXV LT.MINV)G0 TO 20 
DO 16 i=niNv,riHXv 
ANG«U0*i. RANGE'' I J-RANGE'' MINV ) )/C2 
W«CnPLX( cost. f-NG 'i. 3IN( hNG ) ) 
2E«E2+CRO0T(. I ;*FULSEf RANGEC I )-RG >*W 



16 

CONTINUE 



C FILTER 



£0 

X=RHO*r+RH0Ui CABSl 22 ) 

t = X 

IF( IRPTP. EQ» 0 ;C0 T0 1 1 
IPL0T( 1 14QC *X*RLTUD 
rO 24 1= I. IFFTF 



^4 

CALL PLOT IT; I PL 0T , LAST > I ONE- 1 ONE ) 



IS 

FORhHTC IH 4E1S 6- 417 ) 

GO TO 11 

END 



. U 

EXT FUNC 



IPhN 

INT4 VAR 



S 

CMPX VAR 



51) 

CMPX VhR 



5I)D 

CMFX VAR 



CMFLV 

E.7 FUNC 



CONJG 

E’JT FUNC 



W 

CMP>t VAR 



25 

CMPX VAR 



TEL 

REAL VAR 



IGFD 

INT2 VhR 



IPLOT 

INT2 VHP 



LAST 

INT2 VHR 



RhNGE 

FEhL vhR 



CFO'i)T 

PEhL VAR 



THN 

STATE FN 



. Q 

EXT FUriC 



F 

EXT FUNC 



>D 

FORM PHP 



sirj 

EXT FLiriC 



cos 

EnT FUNC 



SPEC 

STATE FN 



'C 

FORM PAR 




REmL VAR 
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EXF 

EXT FUHC 

HhINl 

EXT FUNC 

BGcLI'A 

REAL '<AR 

. R 

EXT FUNC 

FULSE 

STATE FN 

X& 

F0FM PAR 

SIciPLS 

PEAL VAP 

GftIK 

SThTE FN 

XU 

F0RM PAR 

PUG/AL 

REAL VAR 


EXT FUNC 

AES 

EXT FUNC 

SIGTIM 

REAL VAR 

NFFT 

INT2 VAR 

LU 

INT2 VAR 

P2 

PEAL VhR 

C 

REAL VAR 

DT 

REAL VAR 

PULSTh 

REAL VAR 

FO 

REAL VAR 

13 

INT2 VAR 

14 

INT2 VAR 

PTCiR 

REhL VhR 

I0N£ 

INT2 VAP 

lEiG 

INT2 VAR 

IGl'FTS 

1NT2 VAP 

■u‘0 

REAL VAR 

r,‘^FT£ 

INT2 VHP 

NFFTai 

INT2 VAR 


LABEL 

:h 

EXT FUNC 

Zc 

LABEL 

Rflf 

REAL VhR 

Pfl5FHM 

REAL VAR 

PERIOD 

REhL VAR 


LhBEL 

FOTS 

REHL VAP 

31 

LA6EL 

32 

LABEL 

IPFTU' 

INT2 VAF 

IRFTP 

INT2 VhF 

41 

LABEL 

lELKS 

INT2 VAR 

FhO 

REHL VAP 

RhO 1 

PEAL VAR 

SOF^T 

EXT FUNC 

ALT uD 

REAL VAF 

C2 

REAL VAR 

RUGINC 

REAL VhR 

ANGLI.I 

real VHR 

HNGCTR 

REhL VAR 

enLITH 

REAL VAR 
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PNG3DB 

REAL VAR 

fnNGMHA 

REhL VhR 

ANGMIN 

REAL VhR 

to 

REAL VhR 

GDINC 

REAL VAP 

FLOAT 

EXT FUNC 

FLOATS 

EXT FUNC 

5 

LABEL 

I 

INT2 VAP 

. t 

EXT FUNC 

t3I'E 

REAL VAR 

Dt 

REAL VAP 

DtS 

REAL VAR 

t02 

REAL VAR 

DF 

REAL VAR 

RTDF 

REAL VAR 

T>\i 

REAL VAR 

FOUR IT 

EXT FUNC 

IC3UNT 

INT2 VAR 

40 

LABEL 

tF 

EXT FUNC 

1 

LABEL 

A 

REAL VAR 

RANDG 

EXT FUNC 

fM 

EXT FUNC 

I FREQ 

INT2 VAR 

fL 

EXT FUNC 

f A 

EXT FUNC 

SCALE 

REAL VAR 

SCALE 1 

REHL VAR 

SCALES 

REAL VAR 

2 

LAPEL 

X'^ 

REAL VAR 

JK 

EXT FUNC 

3 

LABEL 

IC 

INT2 XAR 

SETGFD 

EXT FUNC 

TMDIF 

PEAL VAR 

REAL 

EXT FUNC 

K 

INT2 VAR 

4 

LABEL 

t 

PEAL VAP 

TANTHT 

REhL VAF 

TkNhLF 

FEhL VAR 

ThNDIF 

REhL VAR 

& 

LABEL 

DLT 

REAL VAR 

DLT 1 

REAL VAF 

HEIGHT 

PEAL ^'HP 

DERIm 

REAL VAR 

DEFIV2 

REAL VAR 

ANG 

REHL VHR 
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AThN 

EXT FUNC 

radcrv 

real VHP 

. A 

EXT FUNC 

30 

LAEEL 

J 

INT2 VAR 

PL0TIT 

EXT FUNC 

NUM 

INT2 VAR 

10 

LABEL 

7 

LABEL 

a 

LABEL 

9 

LABEL 

TIHE 

REAL VAR 


LABEL 

1 1 

LABEL 

RG 

REAL VAP 

RGriH< 

REAL VAR 

rgmin 

REAL VAR 

12 

LABEL 

J3 

LABEL 

. s 

EXT FUNC 

tM\U 

INT2 VAR 

14 

LABEL 

15 

LABEL 

MA\V 

INT2 VAR 

2D 

LABEL 

16 

LABEL 

ChES 

EXT FUNC 

24 

LhBEL 


LABEL 

. 

EXT FUNC 

OOCO EFP0PS 
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Simulated Surface Second Derivative 


Intersection Defines Specular Points* 


Simulated Surface First Derivative 


Simulated Surface 


RANGE > 


OWOWAl, PAoi 


y OF 
PO 



Program SIMUL - Sample Pulse 
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Program 

PUOTIT Do = ‘*jrn:>’ 


PLOTIT 




SCRAT 

ENTRIES AT 3ET* = I i’' ^•-tTIT 

CPLi_-E^E 1^5T CALLED 

FIFS7 T»3 ;e* jp -,FI; F_-:-7I' 5 hLw : = 

CALL SET;f:- KFID' N;,pid rlvS'^ACEJ 

WHERE 

LL»Pl5>>TTER L J 

ItaPID IS AN r-FFHT Cr T iPir LC'-pTnN3 
N^FIf-TpE PEr <_F “/VEIE Frjwr- 
N^'-.PA‘E»NUf*&EP OR F :iTS EFT /EEW 
LIMES IM THE --DIFE'TI^N IF -O CPID LINES 
TO PLOT call 

:ALL PL IDhTF LHS'DHir. NPDIHTS NCONNECT ) 
WHERE IDATfi Ir •- Zfi~ ^ 
rKOtir^EC'"- 0 Tj PL'T jt 5 g^iLr 

N2 T^i C2INWE ^ “«E = r f. t /T'^H 

-HE last if L~ 0> N *E*:' F-”ll O' ?*p TYPE' 
npoirr =E® p^iUTs iii *h7h hi-fbay 

LhSTt-ATP IE PF'^AY " ■ -Zl 2 ’‘-'E l-»ST 
rJlMT* >F I'lrEfJSI N “'Jh. T ripi-IN-: 

L„c7r„-TM 13 .-tL EE n " . imCJ 

IS .E-IFED - r l.f ' 5E E 
hwL ''^L'JEi “'=E ^-EY*I - TZZE^: 

-,.J A'^E H5;-fED 7 t< _'E IE LEE < 0 - 
Pi- P'l- -13-= T-EY _Irl'E* T' 

Thll SRM'SE I'" 'PEY -FE ‘ i EM’^i 


•i^E' 


5Ef '3RD 


entry 

EXTr 1 

STM 

LA 

NHI 


PL F^^TIT 


SETi, 

0 


1 r = F - 
J cT 3 ' 


J _I 


" 

- t 


SI3 

i4 j' 

1 1 

c~ I 


c 


c = 5 

OF 5 



II ~ 


IFF 

L^I 

1 r -- 

1 __ I 

i jFm 

-.I-.O CCOI- 

“ r 


rPL 

S H 

31 

£ 


- 


IC-' 

E 10 

' ; 



*i5 

Zr ^ 

LIS 

4 L 

' ^ cl 

rE 




L.S 

c 

c..*! 

< _ „ 

. ^ r ' 

£ 


L-iI 

'I L I E “ I HE 

c _I 

c, ^ 

-C . 


PI 

w T 

1 . T 1 3 

Jr I 

“I r 

_ 1 ; 



E7 

iC a^ I DC 3 ' 






SIS 


‘ - 

cc J : 


& 


EMP'S 

P’ 


! I 

■ r- 
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LM 

U 4, ,r 

1 i 

Jc 

1 

— w 


LH 

i£ Fi _ J 

^ I 

I'C 

i “ 



STB 

J £ w- 1 1 E i 

' ;i 

Dc 

. M = 

cr 


ITB 

a 'T:=::-k 

JW--I 

•J.EI 

3 1 

£■ V 


LH 

J . 4 

. . -I 

1 ^ - 




EL^S 

4 1 

’ . 7u -I 

c‘Ec 


E Z 

P£ 

El 5 

J C 

* «. 1 I-. I 

c 1 E 


E J 


BM 

P3 

•-EI 

- _' 1 

4e:*i roDo 



LH 

£1 3 14 ' 

n . ^S-I 

cl c: 




Efirs 

F- 

■ rsl 

c4 
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Program RANDG 


tnssn 

RhNI/'i FROG Gf1USSIH^ Ft^NDOH NUMEER GENERATOR LDD Z/76 
EFLST 
JFOFT 

FUNCTION pHNIG(It) 

C It IS hN initialising '/ALUE. LI 9 DIGITS. 

C INTEGEF*4 NOFMhL riUMBEPc GENERATED BY SUMMING 
12 UNIFORMS ilJ The USuhu rASHION, RAI^DU ROUTINE 
C IS USlD not ChlLlI';. 

c 

DO 1 1= 1. 12 
IY=IY*^553S 
IPi. lY. GT . 0 ;Gii 10 2 
IY=IY-^214‘“4S:.5^7+ 1 
2 ^=a+FLOATC IY 

1 CONTINUE 

FH^JDG=X/2 147481.647. - 5. 

Kel uRN 
lND 

RhNDG FuNC '^UE 

RAf^DG FUNC VHF^ 

. 0 E a FUNC 

P E a FUN-' 

IY FORM PhP 

FEHL VAF 
i lhEEl 

1 IfiT4 vhR 

2 Lh&EL 
rLOnT e <T FuNC 


0000 ERP0F5 
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Program RANDU 


fhiSSM 

RhNDU FR&Q random number <aENERAT&»R LDD 3/76 
EFLST 
tFORT 

FUNCTI&N RfiNIU^It) 

C NOTE TriRT IT 15 

C RESULT IS UNIF«?Rh CO> 13 BY lEM 53P 
C lY SHidULD BE INITIALLY ODD V/ITH LT 9 DIGITS 
IY=:Y%655S9 
IFi. lY )5- 6 6 

5 It = IYi-214^4S:6-^7+l 

6 RhND=IY 
FANDU=RHND». 46=c6l3£-5 
RETURN 


PHNDU 

RmNDU 
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lY 
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Rhi'jD 
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END 

FUNC/3UE 
FUNC VAR 
EXT FUNC 
E^T FUNC 
FURM PAF 
LHBEL 
LABEL 
REHL VAR 
EXT FUNC 
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